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Abstract 



Context. Among the components of the infrared and submillimetre sky background, the closest layer is the thermal 
emission of dust particles and minor bodies in the Solar System. This contribution is especially important for current 
and future infrared and submillimetre space instruments - like those of Spitzer, Akari and Herschel - and must be 
characterised by a reliable statistical model. 

Aims. We describe the impact of the thermal emission of main belt asteroids on the 5. ..1000 ^m photometry and source 
counts, for the current and future spaceborne and ground-based instruments, in general, as well as for specific dates 
and sky positions. 

Methods. We used the statistical asteroid model (SAM) to calculate the positions of main belt asteroids down to a 
size of 1 km, and calculated their infrared and submillimetre brightness using the standard thermal model. Fluctuation 
powers, confusion noise values and number counts were derived from the fluxes of individual asteroids. 
Results. We have constructed a large database of infrared and submillimetre fluxes for SAM asteroids with a temporal 
resolution of 5 days, covering the time span January 1, 2000 - December 31, 2012. Asteroid fluctuation powers and 
number counts derived from this database can be obtained for a speciflc observation setup via our public web-interface. 
Conclusions. Current space instruments working in the mid-infrared regime (Akari and Spitzer Space Telescopes) are 
affected by asteroid confusion noise in some specific areas of the sky, while the photometry of space infrared and 
submillimetre instruments in the near future (e.g. Herschel and Planck Space Observatories) will not be affected by 
asteroids. Faint main belt asteroids might also be responsible for most of the zodiacal emission fiuctuations near the 
ecliptic. 
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1. Introduction 

Due to their relatively high apparent brightness at infrared 
wavelength compared e.g. to Galactic stars, asteroids arc 
among the dominant sources at infrared wavelengths and 
can seriously affect the infrared and submillimetre photom- 
etry and source counts. Currently, 377 328 minor planets are 
knowid (as of July 30, 2007) in our Solar System, of which 
about 99% are located in the main belt. On the plane of 
the sky the vast majority of main belt asteroids (MBAs) 
are found at ecliptic latitudes below 20°. They have sizes 
between a few ten meters up to about 1000 km. With tem- 
peratures between 200 and 300 K, the asteroids emit pre- 
dominantly at thermal wavelengths between 5 fim and the 
millimetre range. Deep infrared observations close to the 
ecliptic will ther efore always include some of these mov- 
ing t argets (e.g. iTedesco &: Desertl 120021 : 1 Meadows et al.l . 
|2004| ). Such observations also show that only a small frac- 
tion of the existing minor body population is currently 
known and this population might cause a non-negligible 



Send offprint requests to: Cs. Kiss, pkisscs@konkoly.hu 
^ http:/ /cfa-www. harvard.edu/iau/lists/ArchiveStatistics. html 



confusion noise contribution at certain wavelengths and for 
specific instruments. Tedesco & Desert (2002) measured the 
number of main belt asteroids for the first time in a direct 
way by using the Infrared Space Observatory (ISO). Based 
on a statistical asteroid model the authors concluded that 
there are about 1.2±0.5x 10^ asteroids (> 1 km in diameter) 
in our solar system, twice as many as previously believed. 

A set of celestial sources may affect compact source ob- 
servations in two ways: (1) they contribute to the confusion 
noise, the uncertainty in point source photometry due to the 
fluctuations of the sky background; and (2) they also ap- 
pear as individual sources, that may add 'false' detections 
to source counts. 

Recently, several authors calculated confusion noise 
and detection limits for current/future infrared space 
missions (Spitzer, Akari, Herschel and SPICA). These 
papers considered the two major confusi on noise com- 
ponen ts: the extragalacti c background (jLagache et al.l . 
i2003l: iNegrello et all. I2004D. t he Galactic cirrus emission 
( Kiss et al. . 2005 : Jeong et al.l . l2005f ). or the combination 
of the two (| Jeong et all . 120061 ). 
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It is an important question for infrared space missions 
whether faint asteroids close to or below the detection limit 
could contribute significantly to the confusion noise of these 
instruments or could be present as a significant count of 
contaminating point sources in the field of view. To con- 
sider thes e asteroids, a r e hable statistical model is needed. 
Recently, iTedesco et al.l (|2005[ ) presented the "Statistical 
Asteroid Model" (hereafter SAM). This model is based 
on a population of ^ 1.9 x 10^ asteroids obtained from 
the complete known asteroid sample (as of 1999), plus ex- 
trapolation of the size-frequency distribution (SFD) of 15 
asteroid dynamical families and three background popu- 
lations, to a diameter limit of 1 km. The validity of the 
SAM was demonstrated by com paring SAM pred i ctions 
with ISO measurements at 12 /im iTedesco fc Desertl . l2002[ ) 
and Spi tzer measurements at the 8 and 24 /xm photometic 
bands ([Meadows et al.l . |2004[ ). The asteroid counts from 
both deep surveys show good agreement with the SAM pre- 
dictions. 

In this paper we give (1) estimates for the impact of 
MBAs on infrared (IR) to submillimetre (especially space- 
borne) measurements, such as asteroid number counts and 
confusion noise, for a specific date, sky position, wave- 
length/instrument; and (2) observable predictions for the 
infrared regime which can be used to check the parameters 
of the present model, e.g. size, albedo and positional distri- 
bution. This is the first study that investigates the impact 
of asteroids on infrared and submillimetre photometry and 
source counts in detail, based on a realistic asteroid sample. 

2. Data processing 

2.1. Position calculation 

Orbital elements of all SAM asteroids, including the real 
and predicted ones were calculated for different epochs be- 
tween the time span January 1, 2000 and December 31, 2012 
with an average step size of 5 days. This step size is fine for 
yielding a good coverage of the solar elongation in a year. 
The orbital elements were obtained using accurate numer- 
ical integrations including the effect of all inner and outer 
planets. From the orbital parameters the apparent ecliptic 
coordinates, distances and magnitudes were derived using 
the spatial coordinates of the Earth itself and the abso- 
lute magnitudes (which is known from the SAM database) . 
Some of the spacecrafts, including the Herschel and Planck 
Space Observatories will be located at the L2-point of the 
Earth-Moon system. The difference between the spatial po- 
sition of the barycentre of the Earth-Moon system and the 
prospective position of the L2-point is negligible for the ap- 
parent distribution of the asteroids. However, the Spitzer 
Space Telescope is in a significant distance from Earth, 
therefore we performed the position and all subsequent 
other calculations for the Spitzer Space Telescope coordi- 
nates for a limited time span, covering the expected cryo- 
genic lifetime between January 1, 2004 and December 31, 
2009. The actual positions of the Spitzer Space Telescope 
were taken from the NASA/JPL HORIZONS systenfl. 

To test the positional accuracy of our integration, we 
compared the observable ephemerides yielded by the inte- 
gration and the coordinates returned by the Minor Planet 
& Comet Ephemeris Service of the Harvard-Smithsonian 



Center for Astrophysic^ for a couple of known minor and 
dwarf planets, including Ceres, Pallas, Vesta and Astraea. 
The comparison timespan was almost 105 years, resulted 
by a backward integration from March 6, 2006 to 1901.0. 
The differences between the two sets were always smaller 
than 0?04 for a specific minor planet. We have also com- 
pared this difference for a timespan of 12 years, the same as 
the maximal integration period of the complete simulation. 
In this case the error was definitely smaller than 0?004. 

2.2. Thermal brightness calculation 

For each SAM asteroid and for each date thermal fiuxes 
were assigned at 14 fixed wavelengths {Ao,a:}. These wave- 
lengths were chosen the cover the 5 /xm to 1 mm range 
in a logarithmically equidistant way. If a \i wavelength 
was different from {Ao,fc}, the monochromatic flux values 
were interpolated to the desired for each asteroid, in- 
dividually. For the brightness calcu lations we apphed th e 
Standard Thermal Model (STM, see lLebofskv et al.l . ll986l ). 
In this model the surface temperature distribution is calcu- 
lated using the true heliocentric and geocentric distances. 
The asteroids are described as smooth, spherical and non- 
rotating bodies in instantaneous equilibrium with the so- 
lar radiation. No heat conduction into the surface is con- 
sidered. The correction for beaming, shape and conduc- 
tivity effects is done via the ?7-parameter with a value of 
•q = 0.756. Furthermore, the flux at non-zero solar phase 
angles is obtained by applying an empirical phase correc- 
tion of 0.01 mag deg~^ to the flux calculated at opposition. 
The STM has clear limitations wi th respect to flux accuracy 
(e.g. iMiiller fc Blommaertl . l2004t ) or for modeU ing of minor 
bodies outside the main belt (e.g. iHarri^ . Il998h . but highly 
accurate flux predictions are not crucial for our goals. 

2.3. Number counts, fluctuation power and confusion noise 

Number counts: Two kinds of number count quantities 
are calculated: (1) iVtot(Ai), the total count of asteroids 
in the counting cell, normalized by the solid angle of the 
counting cell (sr~^); and (2) N\i^{Xi, Sum), the count 
of asteroids above the detection limit 5*11,11 in a particular 
counting cell, normalized by the solid angle of the counting 
ceh ilc (sr~i). 

Fluctuation powe rs: The full fluctuation power (see 
iLagache et al.l l2003l for an introduction) is calculated from 
the 'observed' distribution of all asteroids in that specific 
cell, for a specific A^ wavelength: 

In this case the summation runs over all of the asteroids in 
the counting cell (sr^-'^). 

The fluctuation power due to non-detectable asteroids, 
(5i^iini, can be calculated for a specific instrument in a sim- 
ilar way as SFq, but in this case only asteroids below the 
detection limit S'lim are considered: 

6Fnn.{\^,Sun^)= (^) ^ 5|(A,) (2) 



^ http:/ /ssd.jpl.nasa.gov/?horizons 
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Note that Ji^nm is instrument-dependent only through the 
actual sensitivity limit, and other characteristics of the in- 
strument, like spatial resolution, are not taken into account 
in the calculation of ^-Fum- These properties are consid- 
ered in the calculation of the confusion noise. SF^) is a fully 
instrument-independent quantity, and both 5Fq and (5i^iim 
depend strongly on the actual asteroid model. 
Confusion noise: Throughout this paper we assume that 
the local spatial distribution of the asteroids is Poissonian, 
i.e. the same fluctuation power can be used to calculate 
the confusion noise at any spatial frequency, independent 
of the instrument. Thus confusion noise can be calculated 
from the fluctuation powers as: 

(7o(A,,f]p) = (r!p-,5Fo(A,))' ^ (3) 

CriimlAi, S'lim, ^p) = [p.p ■ 5F\ia\[Xi,S\m\)^ (4) 

for the 'full' and sensitivity limit dependent confusion noise 
values, respectively, fip is the effective solid angle of the 
detector, which is not necessarily the physical size of the 
actual pixel/aperture, and the confusion noise applicable 
for detection limits of point sources depends on the point 
sour ce flux ext r action /reconstruction method as well, (see 
e.g. iKiss et ahl . I2005D . "Best estimates" of effective sohd 
angles of various instruments and filters can be found at 
our webpagfl 

Due to their relatively low number density, asteroids in 
our model limit the detectability of point sources through 
the 'photometr i c, rath er than the number density limit (see 
iLagache et al.l . l2003l for a detailed introduction). (Tq and 
aiini are lower limits, since there is an unknown contribution 
of small (fainter) asteroids, which is not considered here (see 
Sect. 4 for a discussion). 

In the following, unless otherwise quoted, we refer to all 
these quantities as specific for a given wavelength or photo- 
metric band and instrument, and therefore the wavelength, 
sensitivity limit and spatial resolution (effective detector 
solid angle) dependences are not marked. 

3. Results 

3.1. General results 

The positions of the asteroids in the original SAM 
model were integrated in the time span January 1, 
2000 to December 31, 2012. At each date a spectral 
energy distribution was assigned to each asteroid us- 
ing the STM, as described in Sect. 12.21 Fluctuation 
powers and number counts have been derived from 
this database, which is pub licly available at the 
URL:" http:/ / kis ag. konkoly. hu/ solarsystem / irsam . A tm j| ' . 
The web-interface Is described m more detail in 
Appendix A. 

For some specific instruments and time spans we con- 
structed SFq, iJi^iim, A'tot, and iViim maps to characterize 
the impact of main belt asteroids for the selected instru- 
ments. An example is shown in Fig. [51 

The fluctuation power and number count maps at dif- 
ferent wavelengths (and the same date) show a very similar 
morphology. Considering the celestial structure, the main 
characteristics are the following (presented in the ecliptic 
coordinate system): 

http:/ /kisag. konkoly.hu/solarsystem/irsam. html 



— The SFq and A^tot distributions (and so SFum and 
Mini) arc symmetric in ecliptic latitude (/3), and show 
a maximum at the ecliptic plane at a specific ecliptic 
longitude (A); SFq, SFum, A^tot and show a similar 
morphology in celestial distribution. 

— The maximum extension in /3 is at the anti-solar point, 
where fluctuation power isocontours form a " bulge" , 
while the minimum is at the celestial position of the 
Sun. When comparing different dates, the bulge around 
the anti-solar point moves along the ecliptic as Earth 
revolves around the Sun, and so does the minimum. 

— If maps presented in the ecliptic coordinate systems (A- 
/3 maps) are transformed to helioccliptic coordinate sys- 
tem ([A — Xq]-P, where A© is the ecliptic longitude of 
the Sun), then the maps of different dates are rather 
similar. The differences are about one order lower than 
the median values of the different maps. Therefore it is 
possible to create time-independent average maps for a 
specific instrument setup, which can give 'average' fluc- 
tuation power and number count estimates, without the 
need to specify the exact day of the observation. 

The main component of the asteroid fluctuation power 
(or confusion noise) as well as the number counts can be 
well represented by a map, which is constant in the he- 
lioecliptic coordinate system. This kind of map can serve 
as a good guideline to characterize the expected impact 
of main belt asteroids for the infrared and submillimetre 
measurements of a specific instrument, if the exact date of 
the observation is not known. We present these maps for a 
handful of instruments in Appendix B, and they are also 
available in FITS format at our webpage^. 

For medium ecliptic latitudes the asteroidal sky changes 
in a timescale of a few weeks, with an actual value depend- 
ing on the instrument (wavelength and sensitivity limit) 
and sky position. Close to the ecliptic plane the amplitude 
of these temporal changes are less pronounced, and there 
are practically no changes at high ecliptic latitudes since 
asteroids are present at these places only sporadically. 

3.2. The importance of main belt asteroids for specific 
instruments 

Asteroids may affect IR and submillimetre observations in 
two ways: they can increase the fluctuation power (and 
hence the confusion noise) level, and contribute to the num- 
ber count of point sources. In our investigated wavelength 
regime the main sources of confusion noise are the extra- 
galactic background and the Galactic cirrus emission. The 
strength of the cirrus emission and confusion noise changes 
rapidly from place to place in the sky and is below the 
extragalatic confusion noise level in the best 'cosmologi- 
cal' windows for most of the space IR instruments (see e.g. 
iKiss et all . [20051) . 

Being constant and present in any direction, extragalac- 
tic background fluctuations represent a minimum value for 
the confusion noise. Therefore we used the respective ex- 
tragalactic fluctuation powers to judge, whether asteroid 
confusion noise has to be considered for a speciflc instru- 
ment. In the calculations presented in Table [1] we consid- 
ered an instrument as " affected" , if the asteroid confusion 
noise level was at least half o f the extragalactic backgrou nd 
component (as calculated bv iLagache et al.Ll2003l 1200^ . 
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In the case of infrared space instruments there is al- 
ways a solar elongation constraint for the actually observ- 
able part of the sky, e.g. 60° < A— A© < 120° for the Herschel 
Space Telescope, 85° < A — A© < 120° for the Spitzer Space 
Telescope and 89° < A — A0<91° for Akari. Instruments 
of these spacecrafts can never look at or close to the anti- 
solar point, where the highest asteroid fluctuation power is 
expected. 

As a general result of the comparison of different wave- 
lengths, instruments working in the far-infrared regime are 
not affected, while instrument working in the 5 /im < A,; < 
30 /xm domain can be severely affected by asteroid confu- 
sion noise, at least along or close to the ecliptic. Although 
far-infrared instruments are not affected by the asteroid 
confusion noise, asteroids above the detection limit may 
have a considerable impact on the source counts. This con- 
tribution can be best estimated for a specific measurement 
(instrument, date, sky position) by the web-interface of our 
infrared asteroid model. 



4. Discussion 

The SAM is limited to asteroids with a lower limit in size 
of 1 km in diameter. There is, certainly, a population of 
asteroids with sizes below this limit. Since the confusion 
noise is calculated using Eq. [21 the impact of very small 
boides (a few hundred meters in diameter and below) is 
minor to the confusion noise. 

To test the effect of small [D < 1 km) asteroids on 
the fiuctuation power and confusion noise we first created 
a size-frequency distribution; for asteroids oi D > 1 km 
we used the number count values estimated by eq. 3 in 
iTedesco et al.l (12005), while for D < 1 km the interpo- 
lated values of iBelton etHI l)1992f ) were used. This SFD is 
a combinat ion of the 'SAM' and 'Galileo team' data points 
in fig. 4 in iTedesco et al.l (|2005f ). A simple relationship of 
S{Xi) oc was assumed between the flux density of an 
asteroid's thermal emission at a specific wavelength S{Xi) 
and its diameter D. Then the specific fluctuation power 
5F{D) was calculated for a _D ± AZ3 interval, using Eq. [T] 
It is clear from Fig. [T] that the fluctuation power is dom- 




log,o(D/l km) 

Figure 1. Relative contribution of asteroids with a speciflc 
size D to the asteroid fluctuation power, using the size dis- 
tribution and the simple flux model as described in Sect. 4. 



{D < \ km) minor planets have a negligible contribution to 
the fluctuation power. Presently orbits for almost 400 000 
asteroids are known, i.e., more than 20% of the 1.9T0^ 
SAM asteroid sample (see Sect. 1), but only about 0.1% 
(~2000 asteroids) have known sizes and albedos which are 
the crucial parameters for the confusion noise estimates. 
The calculated power fluctuations are therefore dominated 
by asteroids which have been extrapolated in the SAM via 
size-frequency distributions. However, knowing only 0.1% 
is far from being statistically relevant and in some parts of 
the sky bright (i.e. large) asteroids are completely missing. 
This implies that confusion noise estimates cannot rely on 
the known asteroid sample alone and the application of the 
SAM model for confusion noise calculations is necessary. 

On the other hand, fluctuation power and confusion 
noise are dominated by the largest/brightest asteroids in 
the actual field, even if the known ones are missing. As 
suggested by Fig. [T] the contribution of Z? < 1 km asteroids 
remains negligible, and a further extension of the SAM to 
smaller diameters would not improve the accuracy of the 
confusion noise calculations significantly. 

Asteroids further out in the Solar System (e.g. trans- 
Neptunian objects) will not contribute considerably to the 
confusion noise or number counts, as deduced from their 
currently known size distribution in Appendix C. 

As discussed in detail in Appendix D, SAM asteroids 
contribute to the absolute brightness of the zodiacal emis- 
sion in a negligible level, however, a set of these asteroids 
may be the dominant sources of zodiacal emission fluctua- 
tion power at some speciflc wavelengths and sky regions. 
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ecliptic longitude (deg) ecliptic longitude (deg) 



Figure 2. Expected distribution of fluctuation power (5Fq, left column) and the total number of asteroids (-/Vtot, right 
column) over the sky of SAM asteroids at 20 ^m on ~ from top to bottom - January 1, March 10 and May 20, 2008 
(ecliptic coordinate system). The position of the solar and anti-solar points are marked with "0" and "O" symbols, 
respectively. 
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Instrument 


A, 


logio(5Fo/Jy^ sr-^) 


A-Aq 


|A/3| 




(^m) 


(EGB) 


(MBA) 


(deg) 


(deg) 


Spitzer/IRAC 


8 


1.49 


2.11 


120 


20 


Akari/IRC 


9 


1.47 


1.21 


90 


13 


Akari/IRC 


18 


1.39 


2.23 


90 


16 


Spitzer/MIPS 


24 


2.04 


2.70 


120 


21 


Akari/FIS 


65 


3.47 


1.34 


90 


- 


Spitzer/MIPS & Herschel/PACS 


70 


3.69 


1.62 


120 


- 


Akair/FIS 


90 


4.06 


0.98 


90 


- 


Herschel/PACS 


100 


4.42 


0.95 


120 


- 


Akari/FIS 


140 


4.42 


0.41 


90 


- 


Akari/FIS 


160 


4.43 


0.23 


90 


- 


Spitzer/MIPS & Herschel/PACS 


160 


4.43 


0.36 


120 




Herschel/SPIRE 


250 


4.29 


-0.21 


120 




Herschel/SPIRE 


350 


3.88 


-0.59 


120 




Herschel/SPIRE 


500 


3.27 


-1.64 


120 





Table 1. Comparison of the extragalactic and asteroid full fluctuation power levels at some representative wavelengths 
of infrared and submillimetre space instruments. Asteroid fluctuation power was derived at the specified X-\q helioe- 
cliptic longitude, and at the ecliptic plane. The columns of the table are: 1) name of the instrument; 2) wavelength; 3) 
extragalactic background fluctuation power (EGB); 4) asteroid fluctuation power (MBA); 5) solar elongation used in the 
MBA fluctuation power calibration, in accordance with the visibility constraints of the instrument; 6) Ecliptic latitude 
range above/below the ecliptic plane, where a specific instrument is considered as 'affected' by asteroid confusion noise. 
Missing values indicate that the instrument is not affected. 



